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ABSTRACT 

We measure the half-light radii of globular clusters (GCs) in 43 galaxies from the ACS Fornax 
Cluster Survey (ACSFCS). We use these data to extend previous work in which the environmental 
dependencies of the half-light radii of GCs in early type galaxies in the ACS Virgo Cluster Survey 
(ACSVCS) were studied, and a corrected mean half-light radius (corrected for the observed environ- 
mental trends) was suggested as a reliable distance indicator. This work both increases the sample 
size for the study of the environmental dependencies, and adds leverage to the study of the cor- 
rected half-light radius as a possible distance indicator (since Fornax lies at a larger distance than the 
Virgo cluster). We study the environmental dependencies of the size of GCs using both a Principal 
Component Analysis ( PCA) as well as 2D scaling relations. We largely confirm the environmental 
dependencies shown in [Jordan et al.l ()2005D . but find evidence that there is a residual correlation in 
the mean half-light radius of GC systems with galaxy magnitude, and subtle differences in the other 
correlations — so there may not be a universal correction for the half-light radii of lower luminosity 
galaxy GC systems. The main factor determining the size of a GC in an early type galaxy is the 
GC color. Red GCs have {rh) = 2.8 ± 0.3 pc, while blue GCs have {rn) = 3.4 ± 0.3 pc. We show 
that for bright early- type galaxies {Mb < —19 mag), the uncorrected mean half-light radius of the 
GC system is by itself an excellent distance indicator (with error ~ 11%), having the potential to 
reach cosmologically interesting distances in the era of high angular resolution adaptive optics on 
large optical telescopes. 

Subject headings: galaxies: ellipticals and lenticular, cD — galaxies: star clusters: general — globular 
clusters: general 



1. INTRODUCTION 

The launch of the Hubble Space Telescope (HST) rev- 
olutionized areas of astronomy which push the limits of 
high resolution imaging. One area which benefited has 
been the study of the half-light radii, r/i, of globular clus- 
ters (GCs). A typical GC has a half-light radius of a 
few parsecs. Ground-based imaging (before the era of 
adaptive optics) could resolve such objects only in galax- 
ies in the Local Group (LG). This limited the statis- 
tics and confined the study to the GC syst ems in dwarf 
and l ate-type galaxies {e.g., the sample in lForbes et al.l 
l2000f l. HST opened the study to GC systems in galax- 
ies out to ~ 30 Mpc, encompassing large numbers of all 
types of galaxies and including both the Virgo and For- 
nax clusters (at 16.5 ± 1.1 Mpc from iMei et all 120071 
and 20.0 ± 1.4 Mpc from iBlakeslee et al.l 120091 respec- 
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tively). 

Interest in the half-light radius of GCs results both 
from the constraints they provide on the formation and 
evolution of GCs and also their possible use as a distance 
indicator, which d a tes b ack to an initial suggestion by 
iShaplev fc Sawveii (fl927l ). The half-hght radii of GCs, 
rather remarkably, a re a lmost independent of GC mass 
(McLaughlin 2000; Jordan et al. 2005; Bar mbv et al. 
2007i: .McLaughlin et al.. .2008.: Harris 2009: Harris et al 
2009bh. at l east to ~ 10^M,t^ plasegan et al.l 120051 : 
Murravl [20091 ). In simulations they have been shown to 
be fai rly constant as the G Cs evolve (Spitzer & Thua3 
1972'; 'Lightman & Shapirol [l97l iMurphv et al.l I1990F 
Aarseth fc Heggie .1998,). and may in fact tr ace the char- 
acter is tic sizes of t he pro to-GC cloud (Mu rrav fc LinI 
[1991 IHarris et al.l I2009bn . In the Milky Way, it 
has long been known that rh increases systematically 
with g alactocentric distance, although HST studies 
(iKundu & Whitmorc 2001; Jordan et alj 120051 : iHarrisI 
have shown that in early types, rh is much closer 
to being constant with galactocentri c radius. 

In this paper we extend the work of lJordan et al.l (|2005l . 
hereafter J05) which used ACS data on GCs in Virgo 
cluster ga laxies (from the ACS Virgo Cluster Survey; 
ACSVCS, Cote ct al.l l2004D to study the environmental 
impacts on the rh of GCs and calibrated a corrected mean 
half-light radius (i.e., corrected for the observed environ- 
mental trends) as a distance indicator. We add to this 
study simi lar data from the A CS Fornax Cluster Survey 
fACSFCS: [Jordan et al.|[2007af) . This adds 43 galaxies to 
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the sample, and also extends the lever arm for calibra- 
tion of the corrected mean half-light radius as a distance 
indicator since the Fornax cluster is at a larger distance 
than the Virgo cluster. 

Other recent works have studied the half-light radius 
of GCs, extending both the total number and range of 
morphologies and environments studied. GC sizes in 
M31 and NGC 5128 were stu died with emphasis o n the 
fundamental plane of GC s bv iBarmbv et al.l (I20Q7D and 
iMcLaughlin et al] (|2008| ). respectively. iBarmby et al.l 
(|2007D show (using data for M31, NGC 5128, the MW, 
the Magellanic Clouds and the Fornax dwarf spheroidal) 
that old GCs appear to have near-universal structural 
properties. Measurements of GC sizes in late type 
galaxies beyond the LG include t he Som brero galaxy 
gpitlcr ct al. 200^ IHarris et al.l l2n09bD. NGC 89 1 
(iHarris et al.l l2n09a'l and NGC 5190 (lForbes"enil[2ni?ih . 
Extremes in host galaxy luminosity ar e encompassed 
by th e work on dwarf galaxies by iGeorgiey et al. 
( 2009 ) a .nd the study of six g iant ellipticals of IHarris 
(pOOai . iDeGraaff et all (|2007l ) studied the SBO galaxy 
NGC 1533 in the Dorado group. At the limit of what 
can be currently done f rom the ground is the work 
IGomez fc WoodlevI (l200l . who used IMACS at the Mag- 
ellan telescopes to study GC sizes in NGC 5128. As we 
will comment below, these works extend many of the re- 
sults on GC sizes we have obtained using the ACSVCS 
and ACSFCS to different host galaxy morphologies and 
environments. All in all, the structural properties of 
GCs seem to share many near-universal properties ac- 
cross galaxy morphology and luminosity, but some dif- 
ferences seem to exist as well. 

The paper is organized as follows. In Section 2 we 
describe the data. In Section 3 we present the distri- 
butions of r/i in Fornax early type galaxies. In Section 
4 we describe the possible environmental dependencies 
which we explore with a principle components analysis 
(PGA) in Appendix A and traditional 2D trends in Sec- 
tion 5. In Section 6 we show the final corrections which 
are discussed in Section 7, both as tracers of the GC for- 
mation and evolution and for implications on the use of a 
corrected mean half-light radius as a distance indicator. 
Section 8 presents a summary of our conclusions. 

2. OBSERVATIONS 

As part of the ACSFCS (| Jordan et al.ll200 7af). a sample 
of 43 early type galaxies in the Fornax cluster were ob- 
served in the F475W (w Sloan g) and F850LP (w Sloan 
z) filters. The ACSFCS sample was constructe d from 
the Forna x Cluster Cata l og (FC C: lFergusonlll989D as de- 
scribed in iJordan et al.l (|2007aD using the FCC galaxy 
morphologies. Half-light radii are measured for all GC 
candidates using the procedure explained in detail in J05 
which fits a PSF-convolved King model to the GCs in 
both bands. The final result quoted here is the aver- 
age of the measurement in both bands. The measure- 
ment of rh is obviously challenging in the cases for which 
the GCs are only marginally resolved. J05 show in their 
Appendix A that their code recovers with no bias 
to r/i ~ 0.3 pixels (which is 0.015" or ~ 1.2 pc at the 
distance of Virgo; ^ 1.4 pc at the distance of Fornax) 
under typical observing conditions, assuming the PSF 
model is correct and that GCs are correctly modeled 
by a King profile. As discussed in J05, the systemat- 



ics undertainties due to the modeling of the PSF are of 
the order of 0.05 WFC pixels (2.5 milli-arcseconds), or 
~ 0.25 pc at the distance of Fornax. The level of sys- 
tematic uncertainties was estimated in J05 by comparing 
the measurements done independently in each of the two 
available bands. This estimate of the level of system- 
atic uncertainty agrees well wi t h ind ependent estimates 
of this quantity given bv lHarrisl(|20Q9l ). who uses different 
datasets and code to measure rh with ACS/ WFC, and 
with those derived from a comparison of measure- 
ments done using two different ACS /WFC data sets of 
M87 using the same code we use here (|Peng et al.i r2009'). 
Thus, the Vh we measure, and in particular their poten- 
tial use as distance indicator, are not tied to our particu- 
lar observational setup beyond possible systematic effects 
at the level of a few milli-arcseconds0 

We follow J05 in the construction of GC catalogs 
in each galaxy. The procedure is described in detail 
in IJordan et al.l ()2009[ ) and summarized in J05, so we 
only list here in brief the cuts made. GCs are se- 
lected using the maximu m likelihood estimated proba- 
bility (| Jordan et al.l 12009 ^ of pcc > 0.5, z-band magni- 
tude, z < 23.35 mag, and colors, 0.6 < [g — z) < 1.7 
mag. This selection on magnitude is roughly equivalent 
to the Virgo cut of z < 22.9 mag at the expected turn- 
over of the GC luminosity function shifted to the greater 
distance of Fornax. Contamination and reliability cuts 
result in a size range of r^, = 0.75-10 pc for objects which 
can be reliably identified as GCs. 

The ACSFCS measures distances using the surface 
brightness fluc tuation (SBF) metho d for all 43 galaxies 
in the sample (iBlakeslee et al.|[2009[ ). therefore allowing 
study of the physical half-light radii of all GCs. We re- 
strict the sample to galaxies with 5 or more GCs (also 
removing FCC 202 whose GC system is overwhelmed by 
its much larger companion, FCC 213). This leaves a final 
sample of 37 Fornax cluster galaxies. Table [1] lists these 
galaxies showing (1) FCC number, (2) B-band magni- 
tude from Ferguson (1989), (3) number of globular clus- 
ters, (4) average value of rh calculated using a biweight 
location estimator, (5) alternative names for the galaxies. 

In some sections, we also consider the sample of GCs in 
Virgo cluster galaxies discussed in J05. Because of slight 
changes in the program which calculated the GC prob- 
abilities, pgCj we find a slightly different sample than 
is used in J05; in general, we find that each galaxy has 
a few more GCs than were used in J05. We also add 
the galaxies VCC21, VCC1833, VCC1440 and VCC1075 
which have very small GC systems and previously were 
just below the cut-off of A^gc = 5 used in J05. 

3. DISTRIBUTION OF Rh IN FORNAX CLUSTER 
GALAXIES 

In this paper we will use several different symbols to de- 
scribe the half-light radius of GCs. We start this section 
by providing a reference table of these symbols. Table [2] 
lists all the symbols used along with a brief description 
of what they represent and a reference to the section of 
the paper that defines it (if appropriate) . 

The distribution of both raw rh, and corrected {fh) 
values for all GCs in the system is shown for the nine 

^ One caveat is that the effects of mass segregation can make the 
measured rjj wavelength dependent, with an expected variation of 
~ 5% between the V- and /-bands (see, e.g..iMadrid ct ah (200.9)). 
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TABLE 1 

Basic Information on Galaxies and GC Systems from 
ACSFCS Sample 



Name B mag A^gc i^h) (arcsec) Other Name 
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most luminous Fornax cluster galaxies in Figure [T] In 
those figures, the solid vertical line shows the biweight 
location estimation of the mean of r^, the dashed line 
is the normal mean and the dotted line is the median 
value. The distribution of in all cases is quite far from 
normal with a long tail to large values of as was pre- 
viously found by JOS for GC systems in the ACSVCS. 
The median value of is therefore always smaller than 
the other estimates of the mean. In addition to galax- 
ies in the Local Group, Virgo, and Fornax, this form of 
the GC size distribution has now been observed in galax- 
ies with a wide range of morphologies and luminosities 
(jGeorgiev et al. 2009; Harris 2009; Harris ct al. 2009b). 

FCC 21 (or Fornax A) shows a rather peculiar distribu- 
tion of GC half-light radii, with a much larger than nor- 
mal number of extended GCs. This skews the average 
for FCC 21 to a much larger value than expected. Note 
that FCC 21 is the most luminous early-type galaxy in 
the Fornax clust er and s hows unmistakable evidence for 
a recent merger (|Schweizer ,1980,; Goudfroo ij 2001, ) which 
may account for this odd distribution of GC sizesEB We 
do not exclude the GC system of FCC21 from the sample 
at this time, but will pay close attention to it in what 

Recently, IDa Costa et al.| 120091) suggested that here may be 
two modes of star formation in dwarf galaxies, a "normal" mode 
with ~ 3 pc and an "extended" mode with ~ 10 pc. This 
"extended" mode, if real, could be potentially related to the excess 
of large GCs we find in FCC 21. 



follows. 

We show in Figure [2] the grouped GC systems of FCC 
21 alone, the remaining bright galaxies {Mb < —19 mag, 
excluding FCC21), and the dimmer galaxies {Mb > —19 
mag), with the fitting function whose parameters are 
given in Equation 27 of J05 overlayed (dashed curve). 
The shape of the distribution for the fainter galaxies 
is very similar to the luminous galaxies, but the mean is 
shifted to a higher value. We also show here the distribu- 
tion of corrected f/j values (red, dot-dashed histogram) 
and for the dimmer galaxies the "extra" corrected, fj^ 
values (blue, dot-dot-dot-dashed histogram) as discussed 
in Section 7. We perform a maximum likelihood fit of 
the function described in Eqns. 24 and 25 of JOS to the 
histograms of for the bright galaxies and fj^ for the 
dimmer galaxies finding best fit values of the parameters 
(/i, /3i, /32 and /) of (2.65, 0.32, 0.18, 0.70) and (2.64, 
0.34, 0.25, 0.61) within a few a of the quoted dispersion 
on these values found from fits to a range of ACSVCS 
galaxies in JOS (whose errors describing the range of fits 
in Virgo cluster galaxies of different colors would there- 
fore be reasonable estimates here too). These functional 
fits are overlayed in Figure[2]as the solid curves. It is clear 
that the fit to the distribution of GC sizes in the bright 
Fornax cluster galaxies is very similar to that found for 
Virgo galaxies by JOS, further demonstrating the uni- 
versality of the distribution of GC sizes in early type 
galaxies, especially the most luminous ones. 

4. FACTORS WHICH MIGHT AFFECT GC SIZE 

It is plausible that many factors may influence the size 
of a GC. This could include properties of the GC them- 
selves (such as mass, metallicity or age) as well as ex- 
ternal factors related both to the local environment at 
the position of the GC, as well as differences in the type 
or size of galaxy in which the GC formed and evolved. 
Following JOS we separate possible dependencies of the 
GC half-light radius into three categories 

1. Internal Factors in which we consider trends of the 
half-light radius with properties of the GCs them- 
selves. We use the z-band absolute magnitude (as 
a proxy for GC mass) and the {g — z) colors which 
correlate with GC metallicity and age. Both quan- 
tities are measured directly from the ACS imaging. 

2. Local Factors in which we consider trends of r^ with 
tracers of the local environment of the GC. We con- 
sider the GCs position in it's host galaxy, using 
the galactic radius scaled to the effective radius of 
the host galaxy (rp/rg) and also the local surface 
brightness, fiz, and local color /ig — fi^ at the posi- 
tion of the GC. The derivation of these quantities 
from A CS images are described in iFerrarese et al] 

iooi. 

3. Global Factors in which we consider trends of r/,, 
with global properties of the host galaxy. For this 
we consider the absolute B-band magn i tudes , AIb 
(from the RC3, Ide Vaucouleurs et al.l (|1991| )). as 
well as the {g — z) color and the average z-band 
surface brightness both measured within the effec- 
tive radius, r^ from the ACS images. The galaxies 
FCC 167 and FCC 26 have very poor fits to their 
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TABLE 2 

Symbols Used to Describe the Half-Light Radii of Globular Clusters. 



Symbol Description 

r'fi Measured half-light radius for an individual GC 

(r^j) Average (usually biweight mean) of all measures of rjj in a given galaxy 

r'f^ corrected for dependence on local galaxy surface brightness (see Section 6.2, Eqn 2) 

(r! ) Average of rj^ 

r'^ r'l^ corrected for dependence of (rj^) on galaxy color (see Section 6.3, Eqns 6 and 7) 

(rj^') Average of rj^' 

r'l^ corrected for dependence on GC color (see Section 6.4) 

(fh) Average of 

f'f^ ffi corrected for a residual dependence of (rj^) on galaxy magnitude (sec Sections 6.3 and 7). 

(f'l^) Average of rj^ 




Fig. 1. — The distribution of uncorrected r^, and corrected ff^ values (in pc) for the nine most luminous galaxies in Fornax. Uncorrected 
values are shown by the forward diagonal hashed histogram, corrected by the backwards diagonal (red) hashed histogram. For the 
uncorrected values, the solid vertical line shows the biweight location estimation of the mean, the dashed line is the normal mean, the 
dotted line is the median value. 



ACS data due to the presence of a massive dust 
disk and a region of star formation the inner re- 
gion of FCC 167 and FCC 26, respectively. In the 
case of FCC 167 no effective radius could be mea- 
sured (so we take the one in the RC3), and in both 
cases the global colors are difficult to define. For 
these galaxies we use the average color from the ar- 
eas of the galaxies used for SBF measurements in 
031akeslce ct al. (2009). 

As is discussed in JOS there is a significant interdepen- 
dency between many of these parameters: for example, 
it is well know that a galaxy's total magnitude correlates 
with color, and the colors of the GCs in a galaxy are 
also correlated with both the galaxy's total magnitude 
and color. This makes the problem an ideal candidate 



for a principal component analysis (PCA) to look for the 
primary correlations. Such an analysis is performed in 
the Appendix A. We will also consider the traditional 2D 
scaling relations below in Section 5. 

5. TRENDS BETWEEN VARIABLES 

In this section we consider separately trends between 
variables, using traditional 2D plots and constructing 
best fit relations. In Appendix A we describe a PCA 
which as discussed in Section 4 above is particularly use- 
ful for learning about the shape of the correlations be- 
tween many different variables when significant intercor- 
relations are expected. However, it is still simpler to in- 
terpret the traditional 2D trends between variables (plot- 
ting one as a function of the other), and by doing this 
successively with different factors the interdependencies 
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Fig. 2. — The distribution of uncorrected (r^ in black solid histogram) and corrected {fh in red, dot-dash histogram) values (in pc) for 
(1) FCC 21, (2) all other bright Fornax galaxies (Mg < — 19mag) and (3) fainter Fornax cluster early types {Mg > —19 mag). Also shown 
for the fainter galaxies is the "extra" corrected, f'^ (in blue dot-dot-dot-dash histogram). The dashed curve shows the function fitted in 
J05 to the distribution of fj, from ACSVCS with its peak normalized to the peak of the observed distribution, while the solid line in the 
two upper panels are our best fit of the same function to the bright Fornax cluster galaxy GC systems, and in the lower panel it is the fit 
to the fainter Fornax cluster galaxy GC systems. 



can also be traced. The PCA in Appendix A has given 
us an idea of the shape of the trends we expect to see 
and which relations we expect to be important and unim- 
portant (showing that the variability in can be well 
described by one factor from each of the internal, local 
and global categories we define in Section 4), however we 
will consider all possible factors again in what follows. 

5.1. Internal Factors 

In this section we consider how properties of the GCs 
themselves might affect the observed sizes of the GCs. 

It has been se en in pre vious HST studies 

(IKimdu fc Whit mord IMH, ,T 05. IBarmbv et al l [2007 : 
McLaughlin et al. 2008; Harris 2009 HHarris et al.ll2009b 



with absolute magnitude. Splitting the GCs in our 
sample into red and blue however we see no significant 
differences in the trend with magnitude. Both sub- 
samples of GCs in early-type Fornax cluster galaxy are 
consistent with having no size-luminosity relation. 

As has also been seen i n previous HST studies of GCs 
in bri ght ellipticals (e.g. iKundu et al.lll99£| |Puzia et ajj 
19991 : iLarsen et al.l[200ll J05. iMadrid et al.l l2009: Harril 



2OO90, we find a clear correlation of rh with GC color, 



Madrid et al that the size of GCs is roughly 



independent of their mass and luminosi ty at least for 
GCs with m asses below M = 2 x 10^ M^ ljHa^egan et all 
120051: iMurravi ^200j). As expected, we confirm this 
property with GCs in Fornax cluster galaxies — we see 
no significant correlation of with the luminosity of 
the GCs (see Figure ^ . We indicate the z-band mag- 
nitude corresponding to M = 2 x 10^ Mp^ (assuming 
M/Lz 1.5 in solar units, iJordan et al.l 12007b'): at 
brighter magnitudes (larger masses) than this there is a 
hint of an upturn in the relation of and magnitude. 
iSpitler et al.l ()2006[ ) suggest that in the Sombrero galaxy, 
an Sa/SO central galaxy in a small group, there is a 
statistically significant trend of the sizes of the red GCs 



with bluer GCs having larger half-light radii (see Figure 
[4]). We follow J05 in restricting the sample in color in 
order to study other correlations, considering blue GCs 
{{g — z) < 1.05) and red GCs {{g — z) > 1.15) separately 
in what follows, and returning to the issue of the color 
dependence later. The reason for separating the GCs by 
color is that many of the properties of GCs are correlated 
with GC color, making it hard to separate the depen- 
dencies on each variable. For example, the dependence 
of (r/i) on GC color will affect the measured dependence 
of (rh) with galactocentric radius through the changing 
ratio of red to blue GCs with this quantity (see §3 in 
J05 for more examples); by restricting the sample by 
color we alleviate these problems. These color cuts cor- 
respond roughly to cuts in metalficity of [Fe/H]< —0.8 
(blue) and [Fe/H]> —0.65 (red), where we have used the 
er npirical rela. tion between [Fe/H] and {g — z) presented 
in iPeng et al.l |2006) derived using Milky Way GC data 
to transform {g — z) color to [Fe/H]. Within each of these 
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Fig. 3. — Half-light radii (r^) vs. z-band absolute magnitude of 
GCs in early type galaxies in the Fornax cluster. Overlayed are 
the median values in bins of 100 GCs. The horizontal line shows 
the median value for the whole sample of 2.94 pc. We indicate the 
z-band magnitude corresponding to Af = 2 X lO*" Mq (assuming 
M/Lz ^ 1.5 in solar units, Jordan ct al. 2007b) where a mass- 
radius relation has previously been shown to emerge. 
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Fig. 4. — Half-light radii (r^) vs. color (g— z) of GCs in early type 
galaxies in the Fornax cluster. Overlayed are the median values in 
bins of 100 GCs. The vertical dashed lines show our color cuts 
for blue GCs {g - z < 1.05) and red GCs {g - z > 1.15). The 
horizontal solid lines show the mean values for the blue and red 
GCs, respectively. 

color bins (shown in Figures |4] and [5]) there is httle ev- 
idence for a trend of with color. This is especially 
true for the red GCs — there may still be a slight trend 
for the bluest of the blue GCs to have large values of 
than the mean. The median value of for blue GCs is 
3.17 pc, and for red GCs it is 2.69 pc, using a bi-weight 
mean we find (r;i)biuc = 3.36 ± 0.03 pc (cr = 1.1 pc) 
and (r^)i.od = 2.83 ± 0.02 pc {a = 1.0 pc). The system- 
atic error on both these measurements at the distance of 
Fornax is ~ 0.25 pc (see JOS, §2). 

Possible explanations for this observed correlation 
between and color are given by projection effects 
combined with a correlation between and galac- 
tocentric distance simi lar to that of the Milky Way 
(jLarsen fc Brodid 120031 ). the combined effects of mass 
segregation and the dependence of stellar lifetimes on 
metallic itv (iJor dan 200^) or different conditions at for- 
mation (|Harri^i2009l) . 

5.2. Local Factors 

Here we will look at how the local environment of a GC 
at the position it is found in its galaxy might affect the 
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Fig. 5. — Distribution of GC colors for the combined sample of 
Fornax cluster galaxy GCs. The vertical solid lines show our color 
cuts for blue GCs {g — z < 1.05) and red GCs {g - z > 1.15). 

observed size of the GC. Using the sample of GCs sepa- 
rated into a blue subsample and a red subsample we now 
look for correlations with projected galactocentric radius, 
Tp, local galaxy surface brightness, fi^ and local galaxy 
color, fig — fiz- Obviously these factors are strongly corre- 
lated. In particular we expect /i oc log[i^(rp)] where F is 
some function describing the shape of the surface bright- 
ness profile of a given galaxy. This does not mean that 
descriptions in the two different variables are equivalent 
since the function F can be quite complicated, however 
there will be significant correlation between the two. In 
ellipticals, it is also generally the case that there is a color 
gradient with a bluer mean color in the outskirts of the 
galaxy than is obs erved in the central re gions (as seen in 
ACSVCS data bv lFerrarese et al.ll2006D . Uhimately we 
will follow J05 in choosing /i^ as the most reliable inde- 
pendent variable, however, to begin with we will consider 
correlations with all three factors. 

In order to place all GC systems on the same scale, we 
normalize their galactocentric radii using the effective 
radii of the galaxy, rg as measured from the ACS images 
(Ferrarese et al 2010, in preparation; note than in FCC 
167 heavy dust obscuration prevented fitting the ACS 
image — in this case is taken from the RC3). Figure 
[HI shows half-light radius, rh vs. projected galactocentric 
distance, rp/r^ for the blue GCs in the nine brightest 
early- type Fornax cluster galaxies. In the two brightest 
ACSFCS galaxies (FCC 21 and FCC 213) GCs are only 
traced to a couple of effective radii, but the vast majority 
of ACSFCS galaxies have GC data out to at least 5re, 
with the maximum value in the sample being lAr^ (for 
FCC 249, re = 7.6"). 

A small trend is visible in most galaxies, with GCs 
at larger galactocentric radii having, on average, slightly 
larger half-light radii. We fit simple linear relations of 
the form log(r;i) = a,- + 6,- log(rp/re) to the data in 
all Fornax cluster galaxies, then take a weighted aver- 
age of these results finding br = 0.029 ± 0.016. There 
is no significant trend of the slope of the relation with 
the extent of the observed GC system {rp^max/fe) sug- 
gesting that we are not hiding a stronger relation by 
being dominated by GCs in bright galaxies for which 
the radial extent of observed GCs is small. In contrast 
to the findings of iSpitler et al.l (|2006f) for the Sombrero 
galaxy, we do not find a significantly stronger trend with 
the red GC population than the blue ones, measuring 



Fornax Cluster GC Sizes 



7 



FCC 21 - FCC 213 



FCC 21! 



NGC 1340 



rCC 1 67 



FCC 276 



^^^^^ 



• FCC 147 



IC 2006 



^ FCC 184 



^ '^itv ^^».- 



012345012345012345 

rp/re Tp/r, fp/r. 



Fig. 6. — GC half-light radius (r;j) vs. galactoccntric distance (in units of the galaxy effective radii) for blue GCs in the nine brightest 
galaxies in the Fornax cluster. The dot-dashed line shows the fits to the individual samples done in log-log space; solid line is the weighted 
mean of the fits to the GC systems in all of the ACSFCS galaxies. These can be compared to a zero trend (dashed line - plotted at the 
mean rh for all blue GCs of rh = 3.2 pc, and the relation which would be seen in the MW GC system after projection effects are considered 
(dotted line). 
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Fig. 7. — GC half-light radius [log(?'h,)] vs. local galaxy surface brightness for blue GCs in the nine brightest galaxies in the Fornax 
cluster. The dot-dashed line shows the fits to the individual samples, solid line is the weighted mean of all the fits to the GC systems in 
all of the ACSFCS galaxies. The dashed line shows a zero trend at the mean of the value of r/j for all blue GCs (3.2 pc). 
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br,rcd = 0.041 ± 0.026, which is the same as the trend of 
the blue GCs to within the la errors. We do the same 
fits for log rh vs /i^ (see Figure [T]) finding in this case 
bfj^ = 0.006 ± 0.003 from a weighted mean of the fits for 
each galaxy. These results are similar to, but a bit flat- 
ter (less significant) than the relations found in JOS for 
GCs in Virgo cluster galaxies of br — 0.07 ± 0.01 and 
bfj, = 0.016 ± 0.003. We do this same fit again for the 
Virgo sample, and find 6^ — 0.011 ± 0.003 (the slight 
difference being due to slight changes in the sample and 
fitting methods). A combination of GC systems in For- 
nax and Virgo galaxies gives b^ = 0.008 ± 0.002. 

We also consider the radii of the blue GC vs. the local 
color of the galaxy {fig — fJ-z), which is shown in Figure 
[8l The outer parts of the early-type galaxies in Fornax 
are all bluer than the inner parts. This corresponds to a 
trend with GC half-light radius such that those with red- 
der local galaxy colors are, on average, slightly smaller. 
However, a fit of the form log(r/i) = Oc + bdpg ~ Hz) 
results in a mean value of be = —0.012 ± 0.013 — a neg- 
ligible trend. 

These results come from the blue GCs only, but no sig- 
nificant differences are seen if just the red GCs are con- 
sidered. The exclusion of GCs in the peculiar GC system 
of FCC 21 (= NGC 1316, Fornax A) which are currently 
included, also has little impact on the final results here. 

We now correct the half-light radii of the GCs to what 
is expected for a GC with an underlying surface bright- 
ness of 21 mag arcsec"^. Since we want a global cor- 
rection for all GCs we will use the slope of trend fit to 
both Fornax and Virgo cluster early-type GC systems 
together of 

r'^ = rhlO-°-^°''<-^'-^'\ (1) 

There is evidence however that the trend is slightly 
stronger in the Virgo galaxies, and slightly weaker (pos- 
sibly zero) in the Fornax galaxies. This is also a weaker 
trend than in J05, who used r',^ = r,^10-"■0^^(^'---'^^) . 

5.3. Global Factors 

We now look for correlations of the GC half-light radius 
with global properties of the host galaxy, such as total B- 
band luminosity and the average surface brightness and 
{g — z) color within the effective radius. We use the half- 
light radius, rj^ corrected for local effects, and consider 
the blue {{g — z)cc < 1.05) and red {{g — z)cc > 1.15) 
GCs separately. 

Figure [5] shows the result for blue GCs in our sample 
of early type Fornax cluster galaxies. All blue GCs are 
shown; the solid circles show (rj^) (biweight mean of r^) 
for each galaxy which has three or more blue GCs. We 
exclude FCC 21 (NGC 1316) from the fits as it is much 
more luminous than all other galaxies and has a pecu- 
liar GC size distribution (as shown in Section 3). The 
dotted line shows the trends with galaxy color and blue 
magnitude observed by J05. We fincF^ 

log(r;,biue.For„ax) =0.515(20) - 0.142(75)[(g - z)gai - 1.5] 
log(r;,biuc.For„ax) =0.500(18) + 0.023(07) (Af^ + 20) 
log(r;,biue,For„ax) =0.534(12) + 0.018(09) ( (a^.) - 19). (2) 

Note that we are using a short-hand notation for errors here 
where 0.515(20) = 0.515 ± 0.020. 



We refit the same relations to Virgo cluster galaxies — 
using a sample that is similar, but not identical, to that 
used in J05 (as discussed in §2) — and find: 

log(r-L,biue,Virgo) =0.453(10) - 0.225(45) [(g - z)gai - 1.5] 

log(r-;,biue,virgo) =0.468(11) + 0.014(04)(Mb + 20) (3) 

Note that as well as this sample being slightly different 
to that used in J05 (possibly of particular importance 
here, some low- luminosity galaxies are added), the cor- 
rection for local effects (a correlation of rh with local 
surface brightness) is also slightly shallower than what 
was used in J05. These difference presumably explain 
the slight changes seen in the relationships we find and 
those reported by J05. 

The slope of the log(r^ j^j^^) vs. galaxy color for For- 
nax galaxy GC systems is very slightly shallower than 
that observed for Virgo cluster galaxies (the same within 
2a); both are consistent with the —0.167 ± 0.054 found 
by J05. We find a trend with magnitude at the ~ 3a 
level in both the Fornax cluster and Virgo cluster galax- 
ies (both the trend and the error on it are twice as large 
in Fornax than in Virgo) , while the slope found for Virgo 
by J05 was consistent with zero within the la errors. We 
attribute this difference to the slightly different samples 
and corrections for local effects as described above. J05 
did not explore trends with the average surface bright- 
ness of the galaxy within the effective radius. In any 
case, we find this trend to be very small — only different 
from zero at about the 2a level. 

Repeating this analysis using red GCs in Fornax/ Virgo 
cluster galaxies, we find slight differences such that the 
trends of size with color and magnitude for red GCs are 
slightly larger than for blue GCs. We still find that the 
trend with color is (slightly) larger in the Virgo cluster 
galaxies, while the trend with magnitude is larger in the 
Fornax cluster galaxies. The trend with average surface 
brightness is still consistent with zero (in fact, the sign 
flips). We also repeat the fits only including GCs within 
5re in the Fornax cluster galaxies (and removing FCC 
213 along which FCC 21 which was already not included) 
to make sure that these trends are not biased by the fact 
that the fainter, bluer galaxies typically have GCs to 
larger radial distances than the more massive ones. This 
cut does not significantly change the fits. 

These differences with the trends of half-light radius 
with global galaxy properties in Fornax early types vs. 
Virgo early types and also with blue and red GCs sug- 
gests that care should be taken applying these correc- 
tions, and further study is needed to determine if there 
truly is a universal correlation. They may also point 
to interesting differences in the formation and evolution 
of GCs in Fornax cluster galaxies versus Virgo cluster 
galaxies and provide more constraints on the differences 
between the blue and red subpopulations of GCs. It is 
possible that different selection procedures for the dwarf 
galaxies included in the ACSVCS versus the ACSFCS 
could cause the differences between Virgo and Fornax. 
However, the impact of selection effects on the galaxies 
on trends of the average sizes of the GCs with global 
galaxy properties is likely to be complicated. 

Figure [TU] shows the same plot as in Figure 9, but for 
blue GCs in both Virgo and Fornax cluster galaxies com- 
bined. The fits in this case are 



Fornax Cluster GC Sizes 



9 



1 .0 
0.8 
0.6 
0.4 

0.2 

0.0 
1 .0 

0.8 

0.6 

0.4 

0.? 

0.0 
1 .0 

0.8 

0.6 

0.4 

0.2 
0.0 



^ FCC 71 


_ • _ 

• FCC 213 

.. .-.1. .;. 


FCC 219 

• -.jr--' 


F 


_ • • 




^ I. NGC 1340_ 


. . FCC 167 _ 

• • • 
• • • • % 


. . FCC 276 


• 




• •. » 




FCC 147 

• •• 
/ • 


IC 2006 

• • 
• • 


FCC "84 

• • 


• ' 


^ i 

•t : 


* y m • • • 

• Vjl^..«» 



1.0 -.2 1.4 1.6 -..8 2.01.0 1.2 1.4 1.6 1.8 2.01.0 1.2 1.4 1.6 1.8 2.0 
/j,^ - |U,2 (mag arcsec"^) /itg - /j,^ (mag arcsec"^) fj.^ — f^^ (mag arcsec"^) 
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Fig. 9. — Top: rj^ vs galaxy color for blue GCs. The filled circles 
show the biweight mean of the distribution for each galaxy which 
has three or more blue GCs. The solid line is fit to these averages. 
The open circle is the outlier FCC21 which is not included in the 
fits. Dotted lines show the fits from JOS to Virgo cluster galax- 
ies. Middle panel: the same, but for rj^ vs galaxy total B-band 
magnitude. Bottom; the same for vs galaxy average surface 
brightness. 



logK,biue,both) =0.473(10) - 0.209(41) [(.9 - z)gai - 1.5] 
log(r-;,,cd, both) =0.409(10) - 0.575(91) [(.9 - z)gai - 1.5] 
log(r-;,biuc,both) =0.480(10) + 0.018(04) (Afs + 20) 
log(r;,,,ed, both) =0.434(14) + 0.039(07)(Mb + 20). (4) 

The addition of the systems in Virgo sHghtly reduces 
the significance of the trend with galaxy magnitude, and 
there is a clear difference in the average values between 
Virgo and Fornax for galaxies fainter than -18 mag (with 
Fornax GC systems having larger average than Virgo) . 
The slope of the trend with galaxy color changes mostly 
because of the offset in zeropoint between Fornax and 



Fig. 10. — As Figure[9]but for blue GCs in both Fornax and Virgo 
cluster galaxies; biweight means for Fornax galaxies are shown as 
circles, for Virgo they are squares. The solid line are the fits to this 
data, dashed lines show the equivalent fits for just the blue GCs in 
Fornax galaxies (as in Figure [9)l. 

Virgo, combined with the distribution of colors of the 
galaxies in the two clusters (i.e., the slope for Virgo is 
parallel to but offset from that in Fornax). 

Following J05 we argue that galaxy color and Mb 
should be almost equivalent as tracers of the global vari- 
ation in r^. We show the color magnitude diagram for 
both ACSVCS and ACSFCS galaxies in Figure HH In 
J05 it was shown that folding the linear variation of log 
with a quadratic relation describing the correlation be- 
tween galaxy color and luminosity in the ACSVCS galax- 
ies reproduces well the mild dependence of on Mb- 
We find indications of a steeper trend of average with 
galaxy magnitude in Fornax, and a shallower trend with 
galaxy color, even though there is no significant differ- 
ence in the color-magnitude diagram. The biggest differ- 
ence in the Virgo and Fornax color-magnitude diagrams 
shown in Figure [Til is that in the Fornax cluster the most 
luminous galaxies are not as red on average as they are 
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Fig. 11. — Color-magnitude diagram for all galaxies in the 
ACSVCS and ACSFCS. Fornax galaxies are shown as filled circles, 
while Virgo are open squares. The relation fit by JOS is shown as 
the dotted line. A quadratic fit to the Fornax data is shown by the 
solid line. The red outlier at Mb ~ —19 mag is FCC184, a face-on 
SBO which may suffer from significant internal reddening. 

in the Virgo cluster (Ferrarese et al. 2006, 2010). It is 
perhaps this slight blueward compression of the colors 
of Fornax early types which means that the color trend 
alone cannot explain all the dependence of GC size on 
galaxy mass. Among the Fornax cluster galaxies, the 
red outlier at Mb ~ —19 mag is FCC184, a face-on 
SBO which may suffer from significant internal redden- 
ing. FCC 167 and FCC 26 (two galaxies with poor fits 
to the ACS images for which we use instead of the aver- 
age color within , the average colo r of re gions used for 
the SBF measurement in IBlakeslee et all f2009l )) are the 
bluest points at Mb ~ —20 and —16.5 respectively. 

In order to fix the correlation of GC sizes with the mass 
of the galaxies, a trend with the galaxy magnitude and 
not just the color is needed, especially in low-luminosity 
Fornax cluster galaxies. However the correlation between 
log rjj and galaxy (g — z) color that we find in both Fornax 
and Virgo cluster galaxy GC systems is consistent with 
J05, so we adopt the same correction to rf[ used there, 

= r;jO°-"[(3-^)-'-i-5l, (5) 

or, combining with the correction for local environment 
(to r',) 

r',[ = r,j0"°-"°^(''^~^^)+°-^'^[(3""^^'5-'"^-^], (6) 

We will return to the issue of a possible extra depen- 
dence of average rj^ value with magnitude in §8. 

5.4. Dependence on GC Color 

We now return to the issue of dependence of GC half- 
light radius on color, [g — z)gc- We use the half-light 
radius corrected for both local and global factors as de- 
scribed above, r^'. 

Figure [T^ shows the corrected half-light radius, rj^' vs. 
GC color for all GCs in all of our sample galaxies com- 
bined. The large squares show averages (biweight means) 
in bins of 0.15 mag of color. The solid line shows are best 
fit to a linear relation of the form log(r^') oc {g — z)qc, 
the best fit slope is & = —0.12 ± 0.01. This is similar 
to the relation found for GCs in Virgo galaxies by J05, 
who found log(r;^') cx -(0.17±0.02)(5 - z)gc) but used a 
slightly different correction for local factors, as discussed 
above. 



Fig. 12. — Corrected half-light radius, r'j^ vs. GC color for all 
GCs in all of our sample galaxies combined. The large squares 
show averages (biweight means) in bins of 0.15 mag of color. The 
solid line shows the best fit linear relation of the form: log{r^'} oc 

(9 - 2)gc- 

J05 contain a lengthy discussion of the correlation of 
GC color with galaxy color and luminosity. Less lumi- 
nous, bluer ellipticals have GC systems which are usually 
metal poor, so have very few red GCs and therefore add 
little leverage to any trends of the GC sizes with GC color 
which will be dominated by the trends in more massive 
galaxies. In order to test that the correlation between 
half-light radius and GC color is the same in all colors 
and luminosities of galaxies we follow J05 in defining sub- 
samples of GCs in galaxies of similar colors. We exclude 
any galaxy which has extremely small numbers of very 
red GCs (they must have a number of GCs with colors 
(5 — z)gc > 1-3 which is equal to at least 10% of the 
blue GCs {g — z)gc < 1-05) since these galaxies add lit- 
tle or no information about the color dependence of the 
half-light radius. We then create bins in galaxy color, by 
requiring at least 80 very red GCs {{g — z)gc > 1-3) in 
each bin; some bins have only a single galaxy. 

The results of this exercise are shown in Figure [131 
All subsamples are consistent with having log(r^') c>c 
— (0.12 — 0.17)(5 — z)gc- Interestingly, there is quite 
a wide range of values of the extent of the galactocen- 
tric radii of the GCs (r^/rg) in these subsets, which 
vary from rp/vg = 0.1-1.2 (in the bin containing FCC 
213) to rp/re = 0.5-8.7 (in the bin containing FCC 63 
and FCC 219). This a rgues against projection effects 
(jLarsen fc Brodij |2003') being the primary mechanism 
explaining the r/i -color correlations. Beyond the large 
ACSVCS and ACSFCS samples, other studies of the de- 
pendence of rh with GC color also suggest that this de- 
pcnde nce is at least partl y intrinsic to the clusters (see 
§3.2 of lHarris et alll2009b l for a s ummary o f this evidence; 
note that an earlier claim by Soit ler et al.l 2006 that pro- 
jection effects are responsible for the r/,-color correlation 
in the Sombrero galaxy has been revised in Harris et al.l 
l2009bf) . Possible mechanisms for creating this i ntrinsic 
depen denc e of size and color have been given in iJordanI 
(|200l andlHarrid (2009). 

6. FINAL CORRELATIONS OF HALF-LIGHT RADIUS 

In this section we list the final correlations of the half- 
light radius of the GCs in Fornax cluster galaxies with in- 
ternal, local and global environmental variables (as listed 
in §4). We find subtly different scaling relations for the 
size of GCs in Fornax cluster galaxies than were found by 
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Fig. 13. — As Figure ri2l but in subsamples separated by galaxy color (as indicated). The dashed hne in each panel shows the fit to the 
entire sample in Figure [12] 



JOS in Virgo cluster galaxies. We find that the trend with 
local environment (measured by local surface brightness) 
may be a bit flatter in the Fornax cluster galaxies than 
in Virgo cluster galaxies (6^, Fornax = 0.006 ± 0.003 as 
opposed to Virgo = 0.011 ± 0.003, while the com- 
bined sample gives — 0.008 ± 0.002). After apply- 
ing this local correction, we look at correlations with 
global properties of the galaxies and find a similar trend 
of average size with galaxy color as was found in the 
GC systems of Virgo cluster galaxies, but there is some 
evidence that an extra correlation with galaxy magni- 
tude might be present in the GC systems of Fornax 
dwarfs. After applying both corrections for local and 
global properties finally we find a slightly shallower trend 
of GC size with GC color (log(r;^'} oc -b.l2(g - z)gc vs 
log(rjj') oc — 0.17(.g — 2;)gc in Virgo cluster galaxy GCs). 
We therefore use corrected radii 

fh = r^lO-0-008(M.-21)+0.17[(g-z)gai-1.5l+0.12[(g-z)GC-1.2] 

(7) 

We also consider an extra correction of 

= raO-°-°2[A^-+2il (8) 

which is only important for galaxies with Mb ^ —19 
mag in the Fornax cluster. 

These subtle differences in the trends raise clear ques- 
tions about the applicability of using the average half- 
light radii of GCs belonging to low-luminosity galaxies 
as a distance indicator. On the other hand, this may 
be hinting at an interesting difference between the for- 
mation/evolution of GC systems in dwarf galaxies in the 
Virgo cluster versus the Fornax cluster — possibly as a 



result of the different environments. 

We provide in Table |3] a listing of the ACSFCS galaxies 
for which we have measurements. Galaxies are ordered 
by absolute B-band magnitude. Also included in this ta- 
ble are the measured color of the galaxy, average surface 
brightness at the position of the GCs, average color of the 
GCs and the average sizes of the GCs (raw, corrected and 
"extra corrected" ) . Figure [HI shows the corrected, {fh), 

extra corrected {r'h), and uncorrected, (r/j) vs. galaxy 
color and magnitude for all galaxies in the Fornax sam- 
ple. The dashed lines show the biweight mean of the 
extra corrected values excluding FCC 21 (dashed line 
show the sigma, also calculated using a biweight location 
estimator). Error bars on the points are the 95% confi- 
dence intervals on the means for each galaxy. It is clear 
that after the corrections are applied the value of (f/j) is 
remarkably consistent over galaxy color and magnitude. 
Nevertheless, two outliers are obvious, having larger than 

usual values of (fh). The first (at Mb 22) is FCC 

21 (NGC 1316 or Fornax A), which was noted earlier as 
a merger remnant showing a peculiar distribution of GC 

half-hght radn. The second (at Mb 19) is FCC193, 

an SO galaxy close to the centre of the Fornax cluster. 
Unlike J05 we do not find that the average rh is consis- 
tent with being constant after just corrections for galaxy 
and GC color and local surface brightness of the disk. 
We find that, in Fornax dwarfs, the average value of 
is larger than expected (as is apparent from the middle 
panel of Figure [H)) . 

7. DISCUSSION 
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Fig. 14. — Corrected, (middle panel; corrections for galaxy and GC color and local disk surface brightness), extra corrected (lower 
panel - including a galaxy magnitude term) and uncorrected, (top panel) vs. galaxy color and magnitude for all galaxies in the sample. 
The dashed lines show the biweight mean of the extra corrected values excluding FCC 21 (dashed line show the sigma, also calculated 
using a biweight location estimator). Error bars on the points are the 95% confidence internals on the means for each galaxy. 



7.1. Half- Light Radius as a Distance Indicator 

Despite the differences in the environmental trends of 
GC half-hght radius examined here, we begin this section 
by re-emphasizing how httle the average sizes of GCs 
actually vary from galaxy to galaxy. The observed trends 
are all very mild and result in only small changes in the 
average half-light radii over significant changes in host 
galaxy color and luminosity. This is particularly true 
in the case of the brightest early-type galaxies in which 
the uncorrected mean half-light radius of the GC system 
already offer a distance indicator with errors comparable 
to those from the best techniques currently available (i.e., 
- 10-15%). 

We now consider in more detail the use of mean half- 
light radius as a distance indicator. We first examine only 
the bright early- type galaxies {Mb < — 19 mag) in which 
we have shown that the uncorrected mean half-light ra- 
dius is very close to constant. Figure [15] shows both 
uncorrected and corrected average rh values for bright 
galaxies {Mb < — 19 mag) in Fornax and Virgo, and also 
NGC 4697 and VCC 575. The values follow very closely 
the expected scaling for a constant mean half-light ra- 
dius. The line is not a fit to the data, but simply this 
relation normalized to a value of {rh) = 0".'033 (solid 
line) at D=16.5 Mpc which is the average values found 
for Virgo giants (both in this work and in J05). 

A simple comparison of the biweight average of all the 
average values of uncorrected rh in GC systems of bright 
Fornax and Virgo early-type galaxies implies a Fornax 



distance of 18.6±1.0 Mpc (based on a Virgo distance of 
16.5 Mpc; Mei et al 2007). Note that the error estimate 
includes only the dispersion in the values of {rh) in both 
Fornax and Virgo systems; the systematic error on this 
from the measurement of the GC sizes is 2.8 Mpc. The 
significant outlier in the Fornax galaxies is FCC 21 which 
as we discussed before has a peculiar distribution of rh 
valueQ Excluding FCC21 from the averaging for For- 
nax we find 19.1±0.9 ± 2.5 Mpc. The dispersion in the 
uncorrected value of {rh) is 0'.'0041 (13%) for all systems 
in Fornax cluster bright early types, but FCC 21 accounts 
for almost half of this value - by excluding FCC21 it drops 
to 0'.'G031 (11%). 

Applying the corrections discussed in §6 and §7 above 
(note we apply our corrections to the GC systems in both 
Fornax and Virgo) we find a mean angular size of 0'.'0293, 
with dispersion 0'.'0041 for all bright Fornax early types 
(0'.'0285 and 0'.'0031 excluding FCC 21) while for Virgo 
we find 0'.'0331 and 0'.'0032. This then implies a Fornax 
distance of = 18.6 ±0.9 ±2.5 Mpc (19.1±0.8±2.7 Mpc 
excluding FCC 21) if the distance to Virgo is 16.5 Mpc. 
The dispersion in the value of the corrected mean half- 
light radius in Fornax early types is 14% (11% excluding 
FCC 21) comparable to the dispersion found by J05 for 
Virgo cluster systems. 

We note that FCC 21 would be readily identified as a pecu- 
liar early-type galaxy in any dataset where the rj, of GCs can be 
measured, as the prominent dust lanes would be clearly seen. An 
rji-derived distance on such a galaxy could be potentially affected 
by a large systematic error. 
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TABLE 3 

Summary of GC sizes in Fornax Cluster, Galaxies 



Galaxy 


Mg 


(9 - Z)Gal 




AT 

"GC 


{(g - z)gc) 


{^h/y ) 


{^h/i ) 




FCC 21 


-22.21 


1.37 


20.07 


1.12 


232 


0.040(02) 


0.038(02) 


0.040(02) 


FCC 213 


-21.00 


1.41 


20.63 


1.26 


698 


0.028(01) 


0.027(01) 


0.027(01) 


FCC 219 


-20.63 


1.51 


21.08 


1.19 


220 


0.029(01) 


0.029(01) 


0.029(01) 


NGC 1340 


-20.33 


1.33 


21.14 


1.00 


137 


0.030(02) 


0.027(02) 


0.026(02) 


FCC 167 


-20.33 


1.09 


20.69 


1.24 


266 


0.032(01) 


0.031(01) 


0.030(01) 


FCC 276 


-19.66 


1.38 


21.51 


1.11 


232 


0.025(01) 


0.023(01) 


0.022(01) 


FCC 147 


-19.56 


1.42 


21.76 


1.15 


201 


0.028(01) 


0.026(01) 


0.025(01) 


IC 2006 


-19.32 


1.44 


21.31 


1.15 


89 


0.034(02) 


0.033(02) 


0.030(02) 


FCC 184 


-19.13 


1.59 


21.44 


1.29 


216 


0.030(01) 


0.031(01) 


0.028(01) 


FCC 83 


-19.12 


1.39 


21.55 


1.08 


173 


0.032(01) 


0.029(01) 


0.027(01) 


FCC 193 


-18.83 


1.39 


21.09 


1.08 


23 


0.046(06) 


0.043(06) 


0.039(05) 


FCC 63 


-18.77 


1.45 


21.91 


1.07 


142 


0.034(02) 


0.032(02) 


0.029(01) 


FCC 170 


-18.71 


1.44 


21.89 


0.98 


30 


0.030(04) 


0.027(03) 


0.025(03) 


FCC 153 


-18.59 


1.39 


21.34 


0.92 


28 


0.037(04) 


0.033(04) 


0.029(03) 


FCC 177 


-18.31 


1.30 


22.82 


0.96 


47 


0.038(03) 


0.032(02) 


0.028(02) 


FCC 249 


-18.20 


1.33 


22.65 


1.00 


102 


0.036(02) 


0.031(02) 


0.027(02) 


FCC 190 


-18.04 


1.38 


22.08 


0.98 


103 


0.041(02) 


0.036(02) 


0.032(02) 


FCC 47 


-18.01 


1.29 


22.62 


1.05 


184 


032(01 "1 


027('01'1 


024f01'l 


FCC 310 


-18.00 


1.32 


20.64 


0.93 


18 


0.028(04) 


0.024(04) 


0.021(03) 


FCC 43 


-17.98 


1.15 


20.81 


1.03 


15 


0.036(06) 


0.030(05) 


0.026(04) 


FCC 148 


-17.90 


1.21 


21.41 


1.02 


13 


0.035(07) 


0.030(05) 


0.026(05) 


FCC 255 


-17.80 


1.22 


22.56 


0.99 


57 


0.035(03) 


0.029(03) 


0.025(02) 


FCC 277 


-17.78 


1.33 


21.61 


1.01 


17 


0.034(03) 


0.030(03) 


0.026(03) 


FCC 55 


-17.70 


1.29 


20.49 


1.08 


12 


0.034(04) 


0.030(04) 


0.026(03) 


FCC 152 


-17.25 


1.15 


21.04 


0.86 


5 


0.050(21) 


0.040(16) 


0.034(14) 


FCC 143 


-17.13 


1.34 


22.83 


0.94 


33 


0.033(04) 


0.028(03) 


0.023(03) 


FCC 95 


-16.88 


1.26 


24.51 


0.99 


13 


0.049(11) 


0.040(09) 


0.033(07) 


FCC 204 


-16.72 


1.23 


22.08 


0.85 


7 


0.050(20) 


0.040(16) 


0.033(13) 


FCC 136 


-16.59 


1.23 


22.62 


0.97 


13 


0.045(05) 


0.038(04) 


0.031(03) 


FCC 182 


-16.56 


1.35 


23.29 


1.02 


22 


0.039(05) 


0.034(04) 


0.028(03) 


FCC 26 


-16.49 


0.64 


22.03 


1.01 


14 


0.041(08) 


0.029(06) 


0.024(05) 


FCC 90 


-16.44 


0.93 


21.71 


0.99 


8 


0.040(10) 


0.030(08) 


0.024(06) 


FCC 106 


-16.39 


1.20 


21.43 


1.07 


5 


0.042(20) 


0.036(17) 


0.029(14) 


FCC 324 


-16.21 


1.17 


22.80 


0.89 


7 


0.036(07) 


0.028(05) 


0.022(04) 


FCC 203 


-16.08 


1.17 


22.76 


0.91 


10 


0.046(07) 


0.036(06) 


0.029(05) 


FCC 100 


-16.07 


1.09 


23.10 


0.94 


9 


0.044(12) 


0.034(09) 


0.027(07) 


FCC 303 


-15.97 


1.13 


22.89 


0.95 


10 


0.051(09) 


0.040(08) 


0.031(06) 



The above discussion relies on an assumed mean dis- 
tance to Virgo to test the reliability of the method. We 
can also use the SBF distances to the galaxies in both 
Virgo and Fornax to provide an updated calibration of 
the method. Using a biweight mean of all the corrected 
mean half-light radii of giant galaxies in both the For- 
nax and Virgo clusters (but excluding FCC 21) we find 
a distance in Mpc of 



D 



0.561 ±0.010 ±0.060 



Mpc. 



(9) 



where (fh) is in arcseconds. This is based on a constant 
value of (fh) = 2.71 ± 0.05 ± 0.25 pc for a GC with color 
{g — z) = 1.2, in a galaxy with color {g — 2:)gai = 1-5 
and at an underlying surface brightness of /iz = 21 mag 
arcsec^^. [Note that one prefers to use uncorrected r^ 
values in giant galaxies in both Virgo and Fornax, the 
result is (r/i) = 2.84 ± 0.07 ± 2.9 pc, so the calibration is 
then D = (0.585 ± 0.010 ± 0.060)/(r?,) Mpc.)] 

Finally, we can also use the geometric calibration 
based on the sizes of GCs in the Milky Way of {fh) — 
2.54 ± O.lpc calculated by J05 to provide a distance 
to the Fornax cluster independent of other extragalac- 
tic distance ladder steps. Using our mean corrected 
(fh) for Virgo and Fornax respectively, we find dis- 
tances of 15.8±0.7 ± 1.4 Mpc and 17.9±1.1 ± 1.8 Mpc 
(18.4±1.0± 1.9 Mpc excluding FCC21). These distances 
are reassuringly close to measurements using other tech- 



niques (see, e.g., the discussion in Blakeslee et al. 2009). 
Note that the quoted errors are statistical, plus the sys- 
tematic error on measuring the half-light radii on ACS 
images. Systematic errors introduced by comparing the 
size of GCs in early types to those in the Milky Way can 
be estimated at about 15-20%. We estimate that the 
best distance to Fornax from this method is 18.4 ± 3.7 
Mpc (where FCC21 has been excluded). 

We also show in Figure [16] the same plot for all galax- 
ies in Virgo and Fornax together (i.e., also including the 
fainter galaxies with smaller numbers of GCs) to illus- 
trate the impact of the environmental corrections on how 
the GC systems in the fainter galaxies fall on the relation. 
If we use all galaxies together we estimate a distance to 
Fornax of I? = 17.2 ± 0.6 ± 1.9 Mpc if no corrections 
are apphed, D = 17.6 ± 0.4 ± 2.3 Mpc if the general 
corrections (for GC color, local surface brightness at the 
position of the GC, and host galaxy color) are applied, 
and D = 20.2 ± 0.6 ± 2.8 Mpc if the extra correction 
on host galaxy luminosity is applied to the sizes of the 
Fornax cluster galaxy GCs. We note that this final mea- 
surement is_v^rY_closeto the 20.0 ± 1.4 Mpc recently 
found bv lBlakeslee et all (|2009D using the SBF method. 

Our prescriptions to correct r^ for it to be useful as 
a standard ruler depend on the particular observational 
setup of our program in that the corrections are given 
using magnitudes measured in g and z. This limits their 
usefulness for other observational setups to some extent, 
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Fig. 15. — Top panel: Average values of (no corrections) vs. inverse distance for the giant galaxies {Mb < —19 mag) in Fornax 
(squares), Virgo (circles). Also includes NGC 4697 (star) and VCC 575 (asterisk). The lower panel shows the same, but for average 
values corrected for dependencies on GC and galaxy color, and local disk surface brightness. The solid line shows the calibration of 
oc 1/D from J05 (normalized to 0'.'033 at D=16.5 Mpc, which is also the average value we find for Virgo giant galaxies). 



but does not preclude their use. First, as we have shown 
above the uncorrected mean half-Hght radius is a good 
distance indicator for bright galaxies {Mb < —190- 
These are the class of galaxies where obtaining a dis- 
tance using rh would be most useful, as they provide 
large number of GCs that are necessary to reduce the 
statistical uncertainty in (r^). In case of a lower luminos- 
ity galaxy observed in other bands, population synthe- 
sis models for early-type galaxies and GCs can be used 
to transform from our bands to the bandpasses under 
consideration. Alternatively, galaxies in well observed 
galaxy clusters such as Virgo, or Milky Way GCs, could 
be used to infer empirical transformations to the chosen 
bandpasses if the appropriate data exist. This is a more 
involved procedure, but as pointed out above, measuring 
distances using (r^) is not ideally suited for low lumi- 
nosity galaxies in any case, due to the small number of 
GCs generally available in such galaxies. Finally, we note 
that our calibration for (rh) hinges on the use of a King 
model to derive the size of GCs, and presumably the use 
of other models could lead to some systematic differences 
in the measured half-light radii. Fortunately, it is possi- 
ble to calibrate away this potential differences by using 
the publicly available data of the ACSVCS galaxies. All 

It can be argued that a distance needs to be known first 
in order to to know if measuring a distance with uncorrected 
is appropriate. But there other ways the nature of a galaxy as 
a luminous elliptical could be established at the required level of 
precision: e.g., by its relative luminosity to other galaxies in the 
case of a cluster of galaxies. 



our measured r^ f or GCs belong i ng to t hose galaxies have 
been published in [Jordan et al.l ([_2009'). If using a differ- 
ent model translates into a systematic difference with the 
rfi measurements presented in [Jordan et alJ ([20091) , this 
difference can be then be taken into account. 

7.2. Half-Light Radii as Tracers of GC Formation and 

Evolution 

It is remarkable how little the mean half-light radius 
of GC systems varies in early type galaxies in both the 
Fornax cluster and the Virgo cluster. The variation is 
only about 20% over half a magnitude change in galaxy 
color and three orders of magnitude in galaxy luminos- 
ity; when only luminous (Mb < — 19 mag) galaxies are 
considered, the variation is smaller still. This property 
clearly provides a significant constraint on models of GC 
formation and evolution. A successful model must also 
account for the red GCs being ~ 20% smaller than blue 
ones, for the lack of dependence of GC size on luminosity 
(or mass), and for the differing trends of GC half-light ra- 
dius with galactocentric distance hinted at in early types 
(where the trend is mild) and late types (where it is much 
steeper). 

In J05, two simple models for the formation of GCs 
were pitted against some of the observational trends of 
the half-light radii. The first was a model in which the 
half-light radius of the GCs is determined by overall pres- 
sure confinement of the proto-GC cloud (e.g., McLaugh- 
lin & Pudritz 1996), while a second model was considered 
in which the radius is a consequence of tidal limitation of 
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Fig. 16. — As Figure [Tsl but now showing both giant and dwarf galaxies, and also including the extra magnitude correction for Fornax 
dwarfs (this is only applied to the GC sizes in Fornax, not in Virgo, NGC 4697 or VCC 575). Dwarfs are shown as open symbols, giants 
as filled. Fornax cluster galaxies are squares, Virgo are circles, NGC 4697 is the star and VCC 575 is the asterisk. We omit error bars for 
clarity — the size of the errors can be seen in Figure 15. 



the GCs. Both models immediately run into a problem of 
predicting a scaling of half-light radius with GC mass (or 
luminosity) which is significantly larger than the almost 
zero trend that is observed (r^ ^ M^^^ and ^ M^/^, 
respectively). Thus, in both cases an unknown mecha- 
nism is required — possibly related to a scaling of star 
formation efficiency with cloud binding energy. Assum- 
ing such a mechanism exists, J05 went on to show that 
in the pressure confined case there is only a mild depen- 
dence of the mean half-light radius with galaxy luminos- 
ity {{rfi) oc L~°-^^), while in the tidal limitation model 
the dependence is exactly zero {i.e. (r/j) ^gai)- They 
therefore argued that both models could predict the neg- 
ligible trend of mean half-light radius with galaxy lumi- 
nosity observed in ACSVCS galaxies. 

In this paper we observe a slight dependence of (r,'j) on 
galaxy luminosity in both Fornax and Virgo cluster early- 
types which clearly differs from the almost zero trend 
observed in the GC systems of Virgo cluster early types 
by J05 (note that slight differences in the Virgo cluster 
GC sample, and treatment of local factors may account 
for this difference). We observe that for blue GCs the 
scaling of size with luminosity is (r^) oc 7j-o o5±o.02 ^^^^ 
for red GCs it is (rj^) oc ]^-0 ii±0.04_ ^ possible explana- 
tion of the difference between the subpopulations could 
be that as the red GCs tend to be concentrated more 
towards the centre of the galaxy they may be more sig- 
nificantly impacted by the size of the galaxy (i.e. the 
depth of its potential well). Overall, this evidence of a 



mild galaxy luminosity dependence of GC sizes, slightly 
favours the pressure confinement scenario over the tidal 
limitation model, especially for the red subpopulation, 
however since the trends are so marginal, and appear to 
depend on the GC selection procedure, we do not con- 
sider this strong evidence. 

We also see some small evidence for a shallower trend 
of the GC sizes with local surface density of the host 
galaxy in the Fornax galaxies than was observed in the 
Virgo galaxies as well as a slightly shallower trend with 
GC color (once global and local correlations are corrected 
for) . Such subtle differences in the environmental depen- 
dencies of the sizes of the GCs in Fornax cluster galaxies 
versus Virgo cluster galaxies may be hinting at a depen- 
dence of the formation mechanisms or more likely the 
structural and dynamical evolution of GCs on the global 
environment of the host galaxies. 

As in J05 the slope of the relation between half-light ra- 
dius and galactocentric distance (br) is significantly lower 
than is found for GCs in the Milky Way (MW), which 
(after projection effects are considered) follow a relation 
with for,MW ^ 0-3 ± 0.15. This study then further con- 
firms the discrepancy between the galactocentric trends 
in early types and in the MW, however we note that 
ou r selection removes GCs with rh > lOpc (as discussed 
in [Jordan et al.l 12009 ^ which could be driving the large 
trend seen in the MW GC system. If real, this differ- 
ence in the GC systems is clearly of interest for further 
study and may point to a difference in the formation or 
dynamical evolutions of GCs in early type galaxies vs. 
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late type galaxies or, possibly, of galaxies in low density 
environments vs. clusters. 

There have been other recent papers which take ad- 
vantage of the ACS on HST to measure the half-light 
radii of GCs in galaxies outside the LG and stud y their 
behav iour with galactocentric distance. .DeGraaff et al.l 
((20071) show that the GCs in NGC 1533 (an SBO galaxy 
in the Dorado group which appears to be transitioning 
from late to early type) have a vs. galactocentric 
radius relation more similar to the GCs in the Milky 
Way than in Virgo or Fornax cluster early type galaxies. 
iSpitler et aP (|2006f ) study the GC system of the Som- 
brero galaxy, an Sa or SO galaxy also in the transition 
region between late and early types and see similar be- 
havior (they find &biue = 0.16 ± 0.04, b,ed = 0.32 ± 0.05 
and foaii = 0.24 zfc 0.03 using the same notation we use 
above). iHarris et al.l ([2009a) studied GC candidates in 
NGC 891 (a nearby edge-on spiral thought to be similar 
to the MW) and while the numbers of GCs in NGC 891 
are quite small the noisy trend seen with galactocentric 
radius seems consistent with that seen in the MW. These 
studies are indicative of a general property of the GC 
systems in late-type galaxies in small groups to have a 
steeper dependence on galactocentric radius than is seen 
in the systems of early-type galaxies in clusters however 
the statistics here are still very small. It would be inter- 
esting to study the GC systems of more late-type galax- 
ies to test this hypothesis. At the moment, it cannot 
be determined if it is the galaxy morphology or, the en- 
vironment in which the galaxy resides, that drives this 
change. A study of GC systems of early types in low 
density environments and late types in higher density 
environments would be instructiv e to disentangle the ef - 
fects of morphology and density. iDeGraafF et al.l (|2007l ) 
tentatively explain the difference as the effect of steeper 
density gradients in small groups like Dorado or the LG 
versus large clusters. Since the sizes of GCs may be lim- 
ited by tidal fields, the argument is that GC sizes should 
have a stronger dependence on radius in small groups 
than in clusters. We find here a (marginally) smaller 
trend of GC radius on galactocentric radius for the GC 
systems in Fornax cluster galaxies than in the Virgo clus- 
ter and while Virgo is the larger cluster, the central den- 
sity of Fornax is about twice tha t of Virgo in term o f 
galaxies per Mpc^ (see Table 1 of iJordan et"aLll2007aD . 
Based on this measure of density we may then be ob- 
serving the trend expected by the iDeGraaff ct al. (2007) 
model. 

It is also interesting that while the size-galactocentric 
radius relation appears to point to differences in the 
GC systems of different Hubble types (or galaxies in 
difference environments), recent studies of the disper- 
sion of the GC luminosity function (GCLF) instead 
favour the idea that GC systems in late type galaxies 
(specifically the Milky Way and M31) are not differ- 
ent from those in early types since they fit on the same 
dispersion-lurninosity relation of systems in early types 
(jJordan et al.ll2007bf ). The color distribution of the GC 
systems of late and early types have also been shown to 
be similar and have similar trends with galaxy luminos- 
ity ([Peng et al.ll2006t) . At the moment the only evidence 
suggesting a significant difference between the GC sys- 
tems in late and early type galaxies is the change in the 
size-galactocentric radius relation. Clearly further work. 



including studies to significantly greater galactocentric 
radii need to be done to confirm this as a general prop- 
erty. 

8. CONCLUSIONS 

We have used data on the half-light radii of GCs be- 
longing to galaxies observed in the ACSFCS to extend 
the work of J05 in which the environmental dependencies 
of the half-light radii of GCs in early-type galaxies in the 
ACSVCS were studied, and a corrected mean half-light 
radius was suggested as a reliable distance indicator. By 
adding data from the ACSFCS, we increase the sample 
size for the study of the environmental dependencies, and 
add leverage to the study of the corrected half-light ra- 
dius as a possible distance indicator (since Fornax is at 
a larger distance than the Virgo cluster). 

We find only subtle differences in the environmental 
dependencies of the sizes of GCs in Fornax cluster galax- 
ies from what was found in the Virgo cluster by J05. 
In addition, we perform a Principal Component Analy- 
sis (PGA) to check that no major correlations are being 
hidden (in Appendix A). 

Looking at 2D relations, we again confirm the well 
known results that there is no correlation between GC 
size, r/i, and mass (for M < 2 x IO^Mq, but that blue 
GCs {{g - z) < 1.05 with (rh) = 3.36 ± 0.03 ± 0.25 pc) 
are about 20% larger than red ones ((g — z) > 1.15 with 
[tk) = 2.83 ±0.02 ±0.25 pc). 

We show that the half-light radii of GCs in early-type 
galaxies in the Fornax cluster increase only slightly with 
galactocentric radius (or decreasing surface brightness) 
as was also found by J05 for systems in Virgo cluster 
early-types. In fact, the trend we find in Fornax clus- 
ter systems is slightly shallower than that seen in the 
Virgo cluster. As was found in J05, the trend of with 
galactocentric distance is significantly shallower (only 2cr 
different from zero in Fornax galaxy GC systems) than 
that observed in our Galaxy perhaps pointing to a differ- 
ent formation scenario, or evolution of GCs in early types 
vs. late types (or in high density regions vs. low density 
regions) . We discuss this observation in light of other re- 
cent studies of the GC systems of late-type galaxies and 
argue there is now some evidence for a general property 
of the GC systems of late- type galaxies in small groups 
to have a stronger trend of half-light radius with galac- 
tocentric distance than is seen in early types in clusters. 
This will provide stronger constraints on the differences 
in formation and evolution of GCs in different types of 
galaxies in different environments. 

We confirm the trend of mean half-light radius with 
galaxy color that was first observed in J05, but show sug- 
gestions that there is a residual correlation with galaxy 
luminosity in the mean half-light radius of GC systems 
of Fornax early-type galaxies which is larger than that 
seen in GC systems of Virgo early types. We revisit the 
two simple pictures of the origin of ru in GCs discussed 
by J05, arguing that the additional trend we observe for 
the mean half-light radius of GC systems in both Fornax 
and Virgo cluster early type galaxies to decrease with 
galaxy luminosity (as (r^) oc i-0-05±o.02 \:,\\x(i GCs 
and (r^) oc 7j-o ii±o.o4 --^^ GCs) may be providing 
some support for the "pressure-confined proto-GC cloud" 
model over tidal limitation of . However since the size 
of this trend appears to depend on the details of the GC 
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TABLE 4 

PCA RESULTS FOR INPUT OF ALL 8 ENVIRONMENTAL FACTORS 



Factor 


VI 


V2 


V3 


V4 


V5 


V6 


V7 


V8 


V9 


cr 




0.15±0.02 


0.18±0.02 


0.69±0.03 


0.17±0.08 


0.40±0.11 


0.47±0.10 


0.13±0.05 


0.07±0.02 


0.02±0.01 


0.157 


{9 - z)gc 


-0.29±0.02 


-0.27±0.02 


-0.29±0.05 


-0.25±0.11 


-0.11±0.18 


0.72±0.11 


0.29±0.07 


-0.16±0.02 


0.03±0.01 


0.231 




0.13±0.02 


0.12±0.02 


-0.49±0.10 


0.68±0.33 


0.18±0.09 


0.26±0.07 


-0.17±0.05 


-O.OliO.Ol 


-O.OlitO.Ol 


0.794 


log(rp/re) 


0.56±0.01 


-0.21±0.03 


0.05±0.02 


-0.07±0.04 


-0.23±0.05 


0.14±0.07 


-0.17±0.03 


-0.15±0.01 


-0.71±0.01 


0.384 


flz Sit Vp 


0.53±0.01 


0.05±0.03 


-0.10±0.04 


-0.28±0.12 


-0.08±0.07 


0.21±0.05 


-0.42±0.07 


0.10±0.02 


0.59±0.01 


1.657 


t^g - fJ-z 


0.25±0.02 


-0.25±0.02 


-0.29±0.06 


-0.26±0.12 


0.75±0.13 


-0.16±0.19 


0.20±0.06 


-0.13±0.01 


-0.04±0.01 


1.413 




0.40±0.01 


0.29±0.02 


-0.17±0.02 


0.00±0.04 


-0.22±0.05 


-0.04±0.08 


0.71±0.11 


0.39±0.02 


0.03±0.01 


1.373 


{9 - 2^)gal 


-0.06±0.03 


-0.64±0.01 


0.06±0.02 


0.12±0.05 


0.03±0.02 


0.01±0.02 


-0.08±0.02 


0.74±0.01 


-O.OOitO.Ol 


0.092 




-0.22±0.03 


0.52±0.02 


-0.22±0.05 


-0.32±0.13 


0.20±0.05 


0.15±0.06 


-0.29±0.05 


0.47±0.01 


-0.38±0.01 


0.927 


E-values (%) 


28.0± 0.3 


21. 6± 0.4 


12. 0± 0.2 


10. 7± 0.2 


9.0± 0.2 


8.2± 0.2 


6.5± 0.2 


3.3± 0.1 


0.7± 0.1 




Ei-\/Ei 




1.3 


1.8 


1.1 


1.2 


1.1 


1.3 


2.0 


4.8 





selection (as it was not present in the similar JOS sample) 
we suggest that it should not be over interpreted. 

We show that for the most luminous galaxies {Mb < 
— 19 mag), the uncorrected mean half-light radius is by 
itself an excellent distance indicator, varying by around 
10-15% across galaxies. This is especially true if we re- 
move the unusual GC system of FCC 21 (Fornax A). 
Once corrected for dependencies on GC and galaxy color 
and local surface brightness we find a constant value 
of (fh) = 2.71 ± 0.05 ± 0.25 pc for a GC with color 
{g — z) = 1.2, in a galaxy with color (g— 2:)gai = 1-5 and at 
an underlying surface brightness of /Lt^ = 21 mag arcsec"^ 
across giant galaxies in both the Virgo and Fornax clus- 
ters (for Virgo alone we find {fh) = 2.67 ± 0.07 ± 0.25 
pc, in Fornax it is (fh) = 2.78 ± 0.12 ± 0.25 pc exclud- 
ing FCC21). The same simple geometric calibration used 
by J05 to estimate an independent distance to the Virgo 
cluster of Dyirgo = 16 ± 2.3 Mpc gives a distance to the 
Fornax cluster of fpornax — 18.4 ± 3.7 Mpc (excluding 
GCs in FCC 21). 

This extension of the work of J05 to include GC sys- 
tems in early-type galaxies in the Fornax cluster adds 
support to the idea of a constant mean in luminous 



early type galaxies, but suggests that the environmen- 
tal dependencies may be subtly different in different en- 
vironments which is especially important in the lower 
mass galaxies. The mean half-light radii in GC systems 
of massive early-type galaxies has the potential to pro- 
vide a geometric distance measurement to bright early- 
type galaxies, which could reach cosmologically interest- 
ing distances (as quantified by J05) in the era of giant 
optical telescopes with adaptive optics. 
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APPENDIX 
PRINCIPAL COMPONENT ANALYSIS 

Principal component analysis (PCA) provides a statistical method to search for correlations between two or more 
correlated variables. Since (as described in §4) the size of a GC may depend on several inter-related properties 
of the host galaxy, local galaxy environment and of the GC itself the problem is an ideal candidate for a PCA. 
PCA decomposes the observed correlations between parameters into a set of eigenvectors and eigenvalues describing 
the main variance seen in the data. This method is now becoming a standard in astronomical applications where 
intercorrelations between several variables are common, so we do not provide a full explanation of the method here, 
but rather refer the reader to standard references below. Here we use the PCA routine in the Astronomv IDL LibrarAjp 
(which performs PCA according to the me thod described in "Multivariate Data Analysis" bv lMurtagh fc Hecklll987t ). 
We follow closely the method described in I Woo et all (|2008[ ) in the interpretation of the PCA results and in the use 
of bootstrap re-sampling to calculate the errors. 

The goals of our PCA will not be to derive expressions for the dependence on on other factors, but rather to 
look at the shape of the trends we expect to see and which relations we expect to be important and unimportant. 
In a sense, we use this to arrive at the optimal number of parameters to describe the rt variations. Traditional 2D 
relations are easier to interpret so we will return to them in the main body of the paper (§5) once we understand the 
overall shape of the expected interdependencies. 

We initially perform a PCA on Fornax cluster GCs using all 8 environmental factors described in Section 4 (in log, 
or magnitude space) as well as logr/j. Numerical results are shown in Tabled In this table (and Tables A2, A3 and 
A4 below) the eigenvector outputs (describing the directions of the main variance in the data) are labelled at V« (for 
a PCA with n variables, i — l..n). We list the projection of each input variable onto these eigenvectors as well as the 
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TABLE 5 

PCA RESULTS FOR INPUT OF REDUCED SET OF 5 ENVIRONMENTAL FACTORS 



Factor 


VI 


V2 


V3 


V4 


V5 


V6 


cr 




0.27±0.03 


0.29±0.04 


0.62±0.04 


0.66±0.03 


0.08±0.05 


0.06±0.02 


0.157 


(9 - 2)gc 


-0.34±0.03 


-0.45±0.03 


-0.26±0.03 


0.64±0.03 


-0.24±0.23 


-0.28±0.02 


0.231 




0.16±0.05 


0.59±0.02 


-0.35±0.03 


0.06±0.06 


-0.54±0.43 


0.12±0.02 


1.657 


Ms - 


-0.25±0.05 


0.49±0.03 


-0.48±0.03 


0.28±0.06 


0.48±0.37 


-0.12±0.02 


1.413 


(9 - ^)gal 


-0.67±0.01 


0.05±0.05 


0.17±0.02 


0.03±0.02 


-0.05±0.05 


0.72±0.01 


0.092 


(Mz> 


0.52±0.02 


-0.33±0.04 


-0.40±0.02 


0.24±0.03 


0.14±0.10 


0.61±0.01 


1.199 


E-values (%) 


29.3±0.5 


24.8±0.4 


17.0±0.4 


12.8±0.3 


10.3±0.3 


5.8±0.2 




E^-l/Ei 




1.2 


1.5 


1.3 


1.2 


1.8 





TABLE 6 

PCA RESULTS FOR INPUT OF 1ST REDUCED SET OF 3 ENVIRONMENTAL 
FACTORS - USING COLORS 



Factor 


VI 


V2 


V3 


V4 


cr 


(g - z)gc 
fig — fiz at rp 

{g - 2)RaI 


0.49±0.02 
-0.56±0.02 
-0.27±0.04 
-0.61±0.01 


0.28±0.06 
-0.40±0.04 
0.84±0.02 
0.22±0.03 


0.81±0.02 
0.37±0.04 
-0.20±0.06 
0.39±0.05 


0.14±0.05 
0.61±0.09 
0.42±0.07 
-0.64±0.10 


0.158 
0.231 
1.657 
0.092 


E-valucs (%) 

E^-l/Ei 


37.0±0.6 


26.3±0.4 
1.4 


20.8±0.4 
1.3 


15.9±0.4 
1.3 





eigenvalue ("E- value") expressed as a percentage of the total. We also show the ratios between successive eigenvalues, 
Ei-i/Ei. The PCA routine performs analysis on standardized variables, meaning that they are rescaled to have a 
zero mean and unit variance. This means that the relative sizes of the eigenvalues tells something physical about the 
relative variance of the data in different "directions". The last column, labelled cr, shows the total variance in the 
input variable before it was rescaled. 

We can see in Table Al that the first two eigenvalues are comparable in size {E1/E2 — 1.3) and about twice as large 
as the third, after which there are 5 vectors of comparable size (with ratios of ~ 1.0 — 1.3), until the last two values 
which are smaller again. In physical terms this can be interpreted as the variance in the data being slightly larger in 
2 directions, then of comparable size until we have two directions which have quite small variance. There there is not 
one much stronger trend between variables than any other. From the size of the components on the first eigenvector 
we see that the strongest trend is between the projected radius (log(rp/re)), the local surface brightness, fj,z (since 
these have the largest components of the first eigenvector). The next strongest trend (two largest projections onto V2) 
is between the galaxy color and average surface brightness. 

Obviously fiz and rp/vg are related through the surface brightness profile of the galaxy, and therefore a relatively 
strong trend between them is expected. Since these trace almost the same physical property (i.e., the local environment 
of the GCs) we chose to drop rp/ve in further analysis (this choice is also discussed in Section 5). 

Running PCA on the remaining 7 environmental factors and the half-light radius shows that the strongest residual 
correlations are between MB,ged and the color of the galaxy. We argue below that M^^gai and {g — z)gai are providing 
complementary information about the host galaxy (because of the well known color-magnitude relation in early-type 
galaxies), and therefore chose to keep galaxy color as a distance independent variable. We will also drop the GC 
magnitude z at this point as the last remaining distance dependent factor. 

Running the PCA on the remaining 5 environmental factors we then find a distribution which shows a relatively 
strong trend between the galaxy color and surface brightness (both local and global), and weaker trends between the 
rest of the variables. The half-light radius {log{rh)) is most strongly projected onto the third eigenvector along with 
(in decreasing order) the color at the position of the GC (/ig — /iz) and the surface brightness at position of the GC 
{^z) — both tracers of the local environment of the GC. The full numeric outcome is shown in Tabled 

We run two final PCA on a set of 3 environmental factors plus the GC half-light radius - one in which the global and 
local environment properties are described by surface brightness, and another using the global and local colors. These 
final PGA on a limited set of parameters now indicate a roughly "spherical" distribution of points, in both of these 
4-D spaces (see numerical results in Tables |6] and [7]) . This indicates that there is no one trend between parameters 
which is significantly stronger than any other. It also seems to not matter if we use the surface brightness or color to 
trace of the global and local properties of the host galaxy. 

The scaling relation between variables X and Y can be found from the primary eigenvector, rescaled by the standard 
deviations of X and Y as 

X oc YipipiY. (Al) 

We do not expect to find exactly the same scaling relations of log here as we will when we consider the 2D trends 
since in the 2D trends we attempt to find a representative fit not simply dominated by the massive galaxies with large 
numbers of GCs (and generally at smaller radii compared to re because of the size of the ACS field). For example 
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TABLE 7 

PCA RESULTS FOR INPUT OF 2ND REDUCED SET OF 3 ENVIRONMENTAL 
FACTORS - USING SURFACE BRIGHTNESS 



Factor VI V2 V3 V4 a 

\ogrh -0.50±0.13 -0.20±0.10 -0.52±0.52 -0.25±0.30 0.158 

(9-z)gc 0.62±0.16 0.06±0.05 -0.05±0.07 -0.50±0.58 0.231 
/i^ at Tp -0.54±0.14 0.03dz0.09 0.46±0.45 -0.34±0.39 1.657 

{fjL^) 0.13ili0.05 -0.96±0.12 0.16±0.13 0.01±0.06 0.927 

E-values (7o) 35.0± 0.5 25.0± 0.2 22.0± 0.4 18.1± 0.4 

1.4 1.1 1.2 

below we often fit the trends in eaeli galaxy separately, then make averages of all the trends in all the galaxies. We 
also consider blue and red GCs separately (although find no significant differences), and consider the trend of logr^ 
with GC color only after the corrections for local and global factors are made. All of this is explained in greater detail 
in Section 5 but is included here to explain why we do not place weight on the scaling relations found from the PCA, 
but rather use it only as a way to make sure we are not missing any unexpected correlations and have included all the 
important factors. 

The final conclusion from our PCA here is that the variability in the values of of GCs can be described reasonably 
well using 3 factors, one each from the internal, local and global properties (as defined in §4). 

REFERENCES 



Aarscth, S. J., & Hcggic, D. C. 1998, MNRAS, 297, 794 
Barmby, P., McLaughlin, D. E., Harris, W. E., Harris, G. L. H., 

& Forbes, D. A. 2007, AJ, 133, 2764 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Blakeslee, J. P., ct al. 2009, ApJ, 694, 556 
Cote, P., et al. 2004, ApJS, 153, 223 

Da Costa, G. S., Grebel, E. K., Jerjen, H., Rejkuba, M., & 

Sharina, M. E. 2009, AJ, 137, 4361 
DeGraaff, R. B., Blakeslee, J. P., Meurer, G. R., & Putman, 

M. E. 2007, ApJ, 671, 1624 
de Vaucoulcurs, G., dc Vaucoulcurs, A., Corwin, H. G., Buta 

R. J., Paturel, G., Fouquc, P. 1991, Third Reference Catalogue 

of Bright Galaxies (Springer: New York) 
Ferguson, H. C. 1989, AJ, 98, 367 
Ferrarese, L., et al. 2006, ApJS, 164, 334 

Forbes, D. A., Masters, K. L., Minniti, D., & Barinby, P. 2000, 
A&A, 358, 471 

Forbes, D. A., Spitler, L. R., Harris, W. E., Bailin, J., Strader, J., 

Brodie, J. P., & Larsen, S. S. 2010, MNRAS, (in press; 

astro-ph/0911:4743) 
Georgiev, I. Y., Puzia, T. H., Hilker, M., & Goudfrooij, P. 2009, 

MNRAS, 392, 879 
Gomez, M., & Woodley, K. A. 2007, ApJ, 670, L105 
Goudfrooij, P., Alonso, M. V., Maraston, C, & Minniti, D. 2001, 

MNRAS, 328, 237 
Harris, W. E. 2009, jVp.I, 699, 254 

Harris, W. E., Mouhcine, M., Rejkuba, M., & Ibata, R. 2009a, 

MNRAS, (in press; astro-ph/0901.3733) 
Harris, W. E., Spitler, L. R., Forbes, D. A., & Bailin J. 2009b, 

MNRAS, 401, 1965 
Ha§egan, M., ct al. 2005, ApJ, 627, 203 
Jordan, A., ct al. 2009, ApJS, 180, 54 
Jordan, A., ct al. 2007a, ApJS, 169, 213 
Jordan, A., ct al. 2007b, ApJS, 171, 101 
Jordan, A., et al. 2005, ApJ, 634, 1002 (J05) 



Jordan, A. 2004, ApJ, 613, L117 

Kundu, A., Whitmorc, B. C. .Sparks, W. B., Macchctto, F. D., 

Zepf, S. E., & Ashman, K. M. 1999, ApJ, 513, 733 
Kundu, A., & Whitmorc, B. C. 2001, AJ, 121, 2950 
Larsen, S. S., Brodic, J. P.,Huchra, J. P., Forbes, D. A., 

Grillmair, C. 2001, AJ, 121, 2974 
Larsen, S. S., & Brodic, J. 2003, ApJ, 593, 340 
Lightman, A. P., & Shapiro, S. L. 1978, Reviews of Modern 
Physics, 50, 437 

Madrid, J. P., Harris, W. E., Blakeslee, J. P., & Gomez, M. 2009, 
ApJ, 705, 237 

McLaughlin, D. E., & Pudritz, R. E. 1996, ApJ, 469, 194 
McLaughlin, D. E. 2000, ApJ, 539, 618 

McLaughlin, D. E., Barmby, P., Harris, W. E., Forbes, D. A., & 

Harris, G. L. H. 2008, MNRAS, 384, 563 
Mei, S., et al. 2007, ApJ, 655, 144 

Murphy, B. W., Colm, H. N., & Hut, P. 1990, MNRAS, 245, 335 

Murray, N. 2009, ApJ, 691, 946 

Murray, S. D., & Lin, D. N. C. 1992, ApJ, 400, 265 

Murtagh & Heck 1987, in Multivariate Data Analysis by 

Murtaghfe Heck [Reidel : Dordrecht 1987], pp. 33-48. 
Peng, E. W., et al. 2006, ApJ, 639, 95 
Peng, E. W., et al. 2009, ApJ, 703, 42 

Puzia, T.H., Kissler-Patig, M., Brodie, J.P., & Huchra, J.P. 1999, 
AJ, 118, 2734 

Shapley, H., & Sawyer, H. B. 1927, Harvard College Observatory 

Bulletin, 852, 22 
Spitler, L. R., Larsen, S. S., Strader, J., Brodie, J. P., Forbes, 

D. A., & Beasley, M. A. 2006, A J, 132, 1593 
Spitzer, L. J., & Thuan, T. X. 1972, ApJ, 175, 31 
Schweizer, F. 1980, ApJ, 237, 303 

Woo, J., Courteau, S., & Dekel, A. 2008, MNRAS, 390, 1453 



